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Abstract 



. We discuss Gaussian fluctuations in a spin glass model with one replica symmetry 

I breaking and we show how non-perturbative fluctuations of the break-point param- 

eter can be included in the longitudinal propagator within linear response theory. 
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■ The aim of this letter is to discuss the fluctuations and more generally the corrections 

to the mean field theory of spin glass models where first order replica symmetry breaking 
occurs. We remind to the reader that replica symmetry can be broken in two different 
ways H]: 

• The function q{x) is discontinuous and it takes only a finite number of values (in most 
cases two). Here the function P{q) is the sum of a finite number of delta functions. 
For example in the case of only one step (IRSB) we have 

(lm{x) = go for X < m, qm{x) = qi for x > m. (1) 

The corresponding function P{q) is given by 

P{q) = m6{q-qo) + {l-m)6{q-qi). (2) 

• The function q{x) is a continuous function and in this case also the function P{q) 
has a continuous part. 
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Some models as the Sherrington Kirpatrick and the Edwards Anderson model belong to 
the second category, other models, the random energy model, Ising spins with p interactions 
and p > 2, the q state Potts model with g > 4, the ROME (random orthogonal matrix 
ensemble) belong to the first category. 

The computation of the fluctuations and the corrections to the saddle point limit is 
rather difficult in the second case, where the form of the propagators is quite involved, and 
requires many powerful tools |@]. 

In the first case (IRSB) the situation was supposed to be much simpler, the propagator 
can be explicitly computed taking care of only the fluctuations of go and qi. The problems 
arise when the fluctuations in the variable m are considered. 

Fluctuations changing m by a small amount are small in some sense and they have to 
be taken into account in the computations, but in some other sense they are large and the 
usual formalism (as we shall see) does not take them into account and must consequently 
be modified. Indeed it is true that when m — m qm{x) — > qm{x) in some sense (for example 
in the norm with finite p) , but the quantity 

sup{qm{x) - qfh{x)) = \qm - qm\oc (3) 

X 

does not go to zero in this limit. 

One of the first results suggesting the necessity of taking care of fluctuations which 
correspond to variations of m is the following [Q. In the random energy model (REM) of 
Derrida the free energy can be written as 

The correct result is obtained as a saddle point in m for large A^. Corrections proportional 
to to the free energy density are clearly connected to fluctuations in m, while if we 
consider the formalism of and we represent the REM as a Ising model with a p-spin 
interaction these corrections cannot come from fluctuations in the q parameters, which 
vanish in this limit. 

More recently it has been shown (see 0) that if one does not take into account the 
fluctuations of m one obtains the wrong result for the specific heat while the correct result 
could be obtained by taking into account the m- fluctuations using a simple (but at this 
stage arbitrary) prescription. 

The aim of this note is to compute part of the fluctuations (the so called longitudinal 
propagator) by using the linear response theory, i.e. by evaluating the variation of the 
function q{x) with respect to an external perturbation. This propagator contains singular 
terms, which are not found using the conventional approach. Correct results for the specific 
heat are obtained using this improved propagator. 

We postpone the computation of the full propagator to a future investigation. Here we 
limit ourselves to the computation of the longitudinal propagator. 
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1 The Model 



We use in our analysis a simple model that we consider representative of the class of models 
with a IRSB saddle point. The model is simply obtained by adding an additional cubic 
term to the usual truncated free energy (W in the following), i.e. 

W[Q] = - hm - f ^Trg^ + ^Trg^ + ^ ^(Q„,)^ + ^ Y^{Q^A , (5) 

where t = Tc — T and Tr stands for trace. We recall that in the SK model a = and 
(3 = 1 (while, for example, in the 3-state Potts model a = 1/2 and (3 is negative 0). In 
the framework of the Parisi ansatz the saddle point of Q is looked for in a subspace in 
which Q can be expressed in terms of a function q{x) defined in the interval [0, 1]. In this 
subspace the functional W[Q] is given by 



\t 1 ( \ a 3 

W\q\ = dx -q\x) - - [xq%x) + 3q^ix) q{y)dyj + -q%x) + j^Q^x) 



(6) 



Below Tc stationarity with respect to the order parameter yields the IRSB solution 

qm{x) = qi9{x - m) (7) 

where the parameters qi and m (go = 0) are obtained by the saddle point conditions as a 
perturbative series in /3, 

r 5 /3r2 25 /JV^ 
?i - + nT. + 



1-a 6(l-a)3 18(l-a; 

/3r 5 
m ~ a H 1 
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1 — a 6 (1 — a 
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The role played by the additional cubic term is to provide a breaking of replica symmetry 
which is located at m ~ a while it is well known that in the SK model m ~ (3t. To 
investigate the stability of this saddle point with respect to Q fluctuations we need the 
eigenvalues of the matrix 
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We find that the eigenvalues of this Hessian which should be positive in order to have a 
stable saddle point, are 



Ao 


— Ai^o 


= — r — gi(m — 1) 


- -/3g?/6 

-^-Pql/Q + qiil-m) 
-f3ql/6 + qim/2 
-l3ql/Q + qi{l-m/2) 

-Pql/Q + qim/2 


Ai 


= Ai,i 


= -T ~ aqi - pqj - qi{m - 2) 




= A2,0 


= — r — gi(m/2 — 1) 




= A2,l 


= -T -aqi- f3qf-qi{m/2-2) - 




= Ao,i,i 


= —T — gi(m — 1) 




= Ai,2,2 


= -r - aqi - (3ql + qi 


Ao,i,2 


= Ao,2,l 


= — r — gi(m/2 — 1) 




= Ao,2,2 


= -T + qi 



(10) 



Using the saddle point values, we find that the minimum eigenvalue belongs to four degen- 
erate sub-families (Aq = Ai^o = Ao,i,i = Ai_2,2) and it is proportional to — /3r^. This shows 
that a coefficient a ^ allows the prescription of keeping a negative coefficient f3 in order 
to have a stable IRSB ansatz (see also [§]), without a negative value for m and without a 
negative eigenvalue. 



2 Fluctuations 

Let us now consider this model in the Gaussian approximation. Our aim is to derive, within 
linear response theory, a longitudinal propagator which takes into account m fluctuations. 
In order to have a consistent check of our computation, at the end of this section we 
shall compare the specific heat obtained through this improved propagator with the usual 
expression obtained through the saddle point solution. 

In the replica approach, the longitudinal propagator can be computed in the discrete 
formulation of replica symmetry breaking, i. e. by using the parameters qo, qi and m and 
the global variations 6qo, 5qi and 5m (as done in p), or by considering the function q{x) 
(see Eq. (^) and the local variations 5q{x) (eventually followed by integration). Clearly, 
if we work in the local formulation, that is the first step toward the analysis in the full 
space, we need a method to deal with the m fluctuations. As previously mentioned, these 
fluctuations induce a non-perturbative (not small in the sense of Eq. (|^)) variations on 
the function q[x) and it is unclear how they can be taken into account in a perturbative 
computation. 

In order to take into account these fluctuations let us introduce in this model an external 
field conjugate with the order parameter: 

W\q] W\q] + C q{x)t{x). (11) 
Jo 

^For a complete and general analysis of the eigenvectors structure in the replica approach, see Q]. In 
what follows we use the notation presented in p|. 
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The perturbation induced by this field on the saddle point solution can be parametrized 
as follows: 

= qi9{x — m — 5m) + 5q{x). (12) 

Within linear response theory, we define the longitudinal propagator by considering the 
response with respect to e, i. e. 

^, . 5qt,{x) ^. , 5m 5q(x) 

The equations for the two components of the propagator follow from the equations of 
motion with the source e 7^ and from their expansion to first order in e, i. e. 



5W[q] 



5q{x 



dW[q,] 



: e{x) (14) 

: -gie(m). (15) 
om 

Using this procedure we manage to consider in a perturbative approach a non-perturbative 
contribution. On the one hand the variations 5q{x) and 5m defined in Eq. (^,0) play a 
different role. The introduction of 5{x — m) as a multiplicative factor of the component 
is crucial in (|T3p because this delta-function separates without ambiguities the two 
contributions. On the other hand the two equations (|1^, |1^) are qualitatively different: 



the first is a functional derivative of the free energy functional W[q] while the latter is a 
derivative of a function W{qe{x)) = Wq^^m{5q{x) , 5m) with respect to 5m. 

There, while in the Eq. (llSD the distribution functions are integrated and no ambiguity 



x—m] 



exists, in the Eq. (0) we have to deal with products of distribution functions (i.e. 9 
6{x — m)5{x — m) and 6'^{x — m)5{x — m)). These products are, at this stage, ill defined 
and a regularization scheme is necessary. In what follows we choose a regularization such 
that 

9''{x — m) = 9{x — m), (16) 
O'^'^i^x — m)5{x — m) = j5{x — m) , (17) 



where the function 5{x — m) that occurs in Eq. ([131) is defined as the derivate of the 
function 6{x — m). Therefore relation (|r^) is the derivative of relation (0), that is the 
only arbitrary choice we make. One can also see that Eq. (p^Tf ) involves the following 
prescription to evaluate the integral of the function q'^{x) on a peaked measure: 

q{x) qi5{x — m)dx = q dq — 



/o Jo k + 1 

By expanding equations ([I^) and (|T5|) to first order in e and by using (|16D and (|l^ we 
obtain following equations for the propagator components: 

-q,e{x - m) I dyj^ + /^i^ + ^q,9ix - m)^ = 5(x - .) 
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1 3 6m 
3 ^^5^ 



—qiS^m — z). 



2'"' Jm 6e{z) 

The corresponding result for the longitudinal propagator ([131) is 

G{x,y) = Go S{x - y) + Gi 6{x - m)e{y - m) + 

- [G^6{x - m) + G^e{x - m)) [G^ 6{y - m) + G^e{y - m)) 

where 



(19) 



(20) 



Gn 



1 



gf/?/6' 



1 



qfP/6 {qlP/6-{l-m)q^] 



3(gf/5/6 



m) qi, 



Gi 



N 



3/4 

\ {qf/3/6 - (1 - m) gi/4)(g2/3/6 - (1 - m) gi) 



^gi (gf/5/6-(l-m) gi/4)' 

(21) 



Two new terms, overlooked by the usual computation, appear in this longitudinal prop- 
agator, the term Gq and G^ . These terms, singular at x ~ m, are the effect of the m 
fluctuations. Let us also note that the asymmetry between the regions x > m and x < m 
in this result is due to the assumption go = on the saddle point. 

To conclude, let us verify the previous result and let us investigate its consequence 
on physical quantities, as the specific heat. It is well known that this quantity can be 
computed through the free energy evaluated at the saddle point or by computing the energy- 
energy fluctuations 0. The computation of a specific heat in the mean field approximation 
through the IRSB saddle point gives 



C(. 



rf2 



W[q] 



SP 



35 



1 — a 



a 



(22) 



where the dependence of the m parameter on the temperature implies a contribution to 
the specific heat also from the variation of m. 

On the other hand, by considering the Gaussian fluctuations at zero-loop order and by 
using our prescriptions to deal with the distribution functions, we find also that 



G{t) = -{I dx dy q\x) q\y)). 



1 ,1 

dx dy qsp{x)qsp{y)G{x,y). 

JO 



(23) 



This shows that the new singular terms, which with our prescription ( [T^ ) produce an 
effect in the computation of the physical quantity (p3|), are necessary to recover the correct 
result. Because of the nature of the x variable in the replica approach, the prescriptions 
for the singular measures are necessary to recover the correct result, while in the discrete 
formalism, where one has to deal with the parameters gi and m only, the regularization is 
not necessary and one naturally recovers the correct result. 

In the case of a continuous breaking of the replica symmetry the longitudinal propa- 
gator computed using the linear response theory do coincide with the one obtained by the 
conventional approach |T^. Our result therefore suggests the following scenario 
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• If we break the replica symmetry in a continuous way by adding an appropriate ex- 
ternal field, the longitudinal propagator is correctly given by the conventional tech- 
niques. 

• If, by removing the external field, the function q{x) becomes discontinuous, the lon- 
gitudinal propagator computed via linear response theory goes to the correct one and 
therefore also the conventionally computed propagator will tend to the same value, 
which is different from the value obtained by applying directly the conventional tech- 
niques. 

• We may only conjecture that a correct computation of all the components of the 
propagator (not only the longitudinal one) may be achieved by using the conventional 
approach after having introduced an external field which breaks the replica symmetry 
in a continuous way and then by sending the external field to zero. 

It is a pleasure to thank Theo Nieuwenhuizen for communicating his results before 
publication and for useful discussion. We also thank David Dean for a careful reading of 
the manuscript. 
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